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Abstract. We report new observations of molecular oxy- 
gen in absorption at z = 0.685 in front of the radio source 
B0218-h357. The hues at 56.3 and 118.7 GHz have been 
observed, redshifted to 33.4 and 70.5 GHz respectively, 
with the 12m at Kitt Peak, 43m at Green Bank telescopes, 
and the 45m Nobeyama radio telescope. Deriving the sur- 
face filling factor of the absorbing dark cloud with other 
lines detected at nearby frequencies, we deduce from the 
upper limits on the O2 lines a relative abundance of molec- 
ular oxygen with respect to carbon monoxyde of O2/CO 
^ 2 10"'^ at 1(7, seven times lower than the previous limit. 
The consequences of this result are discussed. 

Key words: ISM: abundances, general, molecules - 
Quasars: absorption lines - Radio lines: ISM 



1. Introduction 

The interstellar abundance of molecular oxygen has been 
a puzzle for a long time. Since the cosmic abundance of 
oxygen (0/H « 8.5 10"**) is a little more than twice as 
large as that of carbon, at most half of oxygen could be 
locked up in CO, the rest should be in the form of atomic 
01, or in the molecules O2, OH, H2O. The H2O/H2 ratio 
is predicted between 10"*" and 10^^ by chemical models, 
and such predictions are compatible with recent ISO re- 
sults (e.g. Cernicharo et al 1997). Since the abundance of 
OH is only 10^^ or lower, it was long thought that O2 
could be as abundant as CO, i.e. 10~^-10~'* with respect 
to H2 or at least 1/lOth of that of CO (Langer et al 1984, 
Viala 1986). In dark clouds, 01 and O2 are predicted to 
have comparable abundances, while most of the oxygen 
should be atomic in diffuse clouds. To circumvent the at- 
mospheric absorption, the recent searches have used either 
the i^O^^O (Fuente et al 1993, Marechal et al 1997a) or 
the redshifted O2 (Combes et al 1991, Liszt 1992, Combes 



& Wiklind 1995) and have come up with an upper limit 
in the O2/CO ratio of « 1.3 10"^ la. 

We report about further efforts to detect an O2 line 
in absorption, using the redshifted molecular cloud at 
z = 0.68466 absorbing the continuum from the quasar 
B0218-h357 (cf Combes & Wikhnd 1995). The selected 
lines are at lower frequencies, where the continuum flux is 
larger, and should provide a better chance of detection. 

2. Observations 

The observations around 70 GHz have been carried out 
at the NRAO 12-m telescope at Kitt PeakEl in February 
1996. The observations around 30 GHz have been obtained 
with the NRAO 43m telescope in Green Bank in February 
1996, and then at the 45m Nobeyama telescope in Japan 
during several observing runs from February to May 1997. 

Kitt-Peak 12-m Two SIS receivers (68-90 and 90- 
116 GHz) were used at 70 GHz (O2), 87 GHz (CS(3-2)) 
and 105 GHz (HCN(2-1)), with an average system tem- 
perature of 345, 225 and 250 K respectively. The receivers 
were used in dual polarization modes, and tuned in single- 
sideband. The corresponding beam sizes are 98, 79, and 66 
". The observations were done using the wobbler switch- 
ing mode. We used simultaneously broad filterbanks (0.5 
MHz resolution) and narrow band autocorrelators (hybrid 
spectrometer), providing a resolution of 97 KHz. However 
the spectra were smoothed to gain signal-to-noise to veloc- 
ity resolutions of 1-3 km/s. Pointing was checked at least 
every two hours by broadband continuum observations of 
radio-sources; the pointing accuracy was 6 " rms. 

Green-Bank 43— m The receiver was a cooled HFET 



with a system temperature T^ 



120K, and a band- 



width of 500 MHz. The beam size is 56 " at the CS and 
O2 frequencies. The conversion ratio was estimated at 10 
Jy/K. Unfortunately, the efficiency was strongly varying 
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with the hour angle, up to a factor 2-3. The efhciency 
was monitored through observations of strong continuum 
sources (pointing sources) such as 3C 274, 3C 84, and 
the pointing corrections were also obtained on Saturn and 
Jupiter. The observations were made by position switch- 
ing. The autocorrelator backend provided a resolution of 
19.5 KHz, or 0.18 km/s at the O2 frequency. 

Nobeyama 45— m The receiver was a cooled HEMT 
of bandwidth 2GHz and T^^^ 140-170K. The half-power 
beam width of the telescope is 59 " and 51 " at 29 and 
33 GHz respectively. The aperture efficiency is 60% at 
30 GHz, leading to a conversion ratio of 2.9 Jy/K, for 
the antenna temperature T^. The observing procedure 
was position switching. The was checked every one to two 
hours by observing a nearby SiO maser star, S per, and 
its accuracy was better than 10 " (peak value). The 2048- 
chanels AOS backends provide a resolution of 37 kHz, i.e. 
0.38km/s at the CS frequency and 0.33 GHz at the O2 one. 
The flux density of the continuum emission was measured 
with a continuum backend in beam-switching mode at 29 
GHz and with AOS at 33 GHz simultaneously with the 
O2 hne. 




-5C 

Velocity (km/s) 



Fig. 1. From top to bottom: HCN(2-1) spectrum at IRAM 
July 1996, HCN(2-1) at Kitt-Peak February 1996; CS(l-O) at 
Nobeyama February-May 1997, CS(l-O) at Green-Bank Febru- 
ary 1996. Spectra have been normalised to the continuum seen 
by the molecular cloud, i.e. 33% of the total continuum. 



3. Results 

The BL Lac object B0218-I-357 is remarkable by its ra- 
dio images, that reveal two distinct flat-spectrum cores 
(A and B component), with similar spectrum and polar- 
ization level, strongly suggesting a gravitational lens phe- 
nomenon (Patnaik et al 1995). That A and B are gravita- 
tional images is confirmed by the presence of a small Ein- 
stein ring surrounding the B image, of 335 milli-arcsecond 
in diameter (Patnaik et al 1993). Since the ring has a 
steeper spectrum, it is probably the image of a hot spot 
in the jet component, and its contribution in the millimet- 
ric domain should be negligible. The Icnsing galaxy at a 
redshift z — 0.68466 absorbs the radio continuum of the 
remote quasar in the HI line (Carilli et al 1993), and in 
molecular lines (Wiklind & Combes 1995). The intensity 
ratio between the two images is A/B k. 3-4 at several radio 
wavelengths (O'Dea et al. 1992, Patnaik et al. 1993). As in 
the similar absorbing lens galaxy in front of PKS1830-211 
(cf Wiklind & Combes 1996), some of the main molecu- 
lar lines are optically thick as attested by the detection 
of isotopomers, but the radio continuum is not totally ab- 
sorbed, which means that the absorbing material does not 
cover the whole surface of the continuum source. In the 
case of PKS1830-211, interferometric observations have 
confirmed that only one image of the source is covered 
by the absorbing cloud (Wiklind & Combes 1997a). In 
B0218-I-357, the situation is still unclear, although it has 
been claimed that only part of the A image is covered at 
cm wavelengths (Menten & Reid 1996) and through op- 
tical imaging (Grundahl & Hjorth 1995). At the period 
of our observations, the fraction of the total continuum 
absorbed was determined to be one third. 

To derive the optical depth of a given line, or at least an 
upper limit on the optical depth, we must know exactly the 
total continuum flux of the source at the epoch of obser- 
vations, since it might be variable, and the surface filling 
factor of the absorption, i.e. the fraction of the continuum 
source surface covered by molecular clouds. The best is 
even to bypass these two steps, by observing another line 
at a neighbouring frequency, which is optically thick and 
gives directly the absorbed continuum level. This was done 
in our previous work, where the neighbouring CO transi- 
tion was shown to be highly optically thick through ob- 
servations of i^CO and C^^O (Combes & Wiklind 1995). 

To calibrate the N,J = 1,1 - 1,0 118 GHz O2 line, 
redshifted at 70 GHz, we could not observe the 115 GHz 
CO(l-O) line unfortunately, since the 12~m receiver was 
not operating at 68 Ghz. Instead, the line CS(3-2) and 
HCN(2-1) were observed at redshifted frequencies of 87 
and 105 GHz respectively. The HCN(2-1) profile is dis- 
played in figure |l|, normalized to the absorbed continuum. 
It can be compared with the analogous spectrum obtained 
with the IRAM telescope in July 1996, also normalised 
to the continuum. Both appear optically thick. From the 
millimetric observations carried out with the IRAM-30- 
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Fig. 2. Radio continuum spectrum of B0218+357. Measure- 
ments from the literature (see NED for a compilation) are plot- 
ted with filled triangles, and the stars correspond to the present 
work. The dashed line is not a fit, but represents a power law 
of slope —0.25. 



Fig. 3. Negative results found for the two selected O2 lines: 
Top : the N,J = 1,1 - 1,2 line observed at Nobeyama at 33 
GHz; Bottom : the N,J = 1,1 - 1,0 line observed at Kitt Peak 
at 70 GHz. Spectra have been normalised to the continuum 
seen by the molecular cloud, i.e. 33% of the total continuum. 
The velocity resolution is smoothed for both spectra to 3.6 
km/s. 



m telescope in 1994, 1995 and 1996, the total continuum 
does not seem to have varied, at least within the error bars 
(10-20%). The CS(3-2) was not detected with the 12-m 
(although it is detected with the IRAM-30-m, Wiklind & 
Combes 1997b). 

To calibrate the N,J = 1,1 - 1,2 56 GHz O2 hne, red- 
shifted at 33 GHz, we have observed the CS(1~0) line, red- 
shifted at 29 Ghz, together with the corresponding ^'^CS 
Unc. The CS(l-O) Hne was detected both at NRO and GB 
(see figure ||), but are not saturated, and the ^^CS line 
was not detected, with an rms noise level of 4mK in chan- 
nels of 4km/s (see Table 1). We must therefore rely on 
the measurement of the total continuum, and on the as- 
sumption that the fraction of the surface covered by the 
molecular cloud is the same at 1cm and 3mm. 

The total continuum level measured at the various 
frequencies and telescopes is displayed in figure ||, to- 
gether with some other measurements in the literature 
(see Combes & Wiklind 1997 for details). Within the er- 
ror bars, aU measurements follow a power-law distribution 
as a function of frequency. In February 1996, the total con- 
tinuum level measured at 30 GHz at Green-Bank was T^ 
= 70 ± lOmK (0.7 Jy), and in February 1997 it was 0.255 
± .017 K at 33 GHz and 0.245 ± .020 K at 29 GHz at 
Nobeyama (i.e. 0.74 and 0.71 Jy) in very good agreement. 

We have assumed thermal equilibrium; this is justi- 
fied since the absorbing molecular cloud is dense with a 
large optical depth, and the O2 lower levels are completely 
thermalised in these conditions (see Marechal et al 1997b). 
Then, to estimate column densities, the usual formula can 



be used (cf Combes & Wiklind 1995), assuming gaussian 
line profiles, with the intrinsic line-width of 4.6 km/s. 



Table 1. Observed molecules 



Mol. 


Transition 


^rcst 


"obs 


Dot.") 


Tel.=) 






GHz 


GHz 






CS 


0^1 


48.990978 


29.080632 


Yes 


NRO 


CS 


0^1 


48.990978 


29.080632 


Yes 


GB 


CS 


2^3 


146.969033 


87.239579 


No 


KP 


l^CS 


0^1 


46.247567 


27.452167 


No 


NRO 


HON 


1^2 


177.261111 


105.220700 


Yes 


KP 


O2 


lo ^ ll 


118.7503 


70.489179 


No 


KP 


O2 


I2 ^ ll 


56.26476 


33.398288 


No 


GB 


O2 


I2 ^ ll 


56.26476 


33.398288 


No 


NRO 



a) Rest— frequencies from the JPL eatalog (Pickett ct al) 

b) Detection: Yes/No. 

c) Telescope: KP=Kitt Peak 12-m, GB= Green-Bank 43-m, 
NRO=Nobeyama 45-m. 



From the observed profiles of CS(l-O) and CS(3-2), a 
column density of N(CS) of 1.5 lO^'' cm~^ is found, with 
an excitation temperature of T^; = 6j^2 K, very close to the 
background temperature at the considered redshift Tbgo 
(l-fz) — 4.59K. With the previously derived H2 column 
density of 5 10^3 cm-^ (Combes & Wikfind 1995), the 
CS abundance is estimated to 0.3 10~^. With this low 
abundance, no signal could be expected for the ^'^CS line. 

The O2 profiles at 33 and 70GHz are shown in figure |[ 
At 33GHz, the NRO upper limit is more significant than 
the Green-Bank one, and we show only that one. They 
correspond to an upper limit of optical depth at Ict of r ^ 
0.025 and 0.12 respectively. According to the adopted ex- 
citation temperature, upper limits on the O2 abundance 
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are derived in Table 2. They arc compared with the col- 
umn densities derived from previous observations of the 
CO molecule (IRAM April 1995), assuming unsignificant 
time variation (Combes & Wiklind 1995). 

The present data together with those obtained earlier 
(Combes & Wiklind 1995) cover a wide range of conditions 
in the absorbing gas, specifically low and high excitation 
temperatures. It is therefore not possible to have a large 
abundance of O2 'hidden' at certain rotational levels. This 
is an important aspect of measuring the O2 abundance 
since the molecular structure inhibits coUisional excitation 
of the N,J = 1,0 level from the N,J = 0,0 level (the 118.7 
GHz line). This docs not apply to radiative transitions 
(Bergman 1995). However, we can see in Table 2 that the 
O2/CO abundance ratio does not depend much on the 
adopted temperature. We find for any choice of that 
the O2/CO ratio is < 2. lO'^ at la. 

4. Discussion 

The new O2 limit derived above on the abundance ra- 
tio O2/CO < 2 10~^ begins to constrain significantly the 
chemical models. Marechal et al (1997b) have recently re- 
considered the chemistry and rotational excitation of the 
O2 molecules in interstellar clouds, cither dark, translus- 
cent or diffuse. They have used the most recent determi- 
nations of reaction rates and collisional cross-sections, and 
explored the effect of UV radiation, density and temper- 
ature. They consider only steady-state equilibrium. Their 
conclusions show that the most important parameters de- 
termining the O2 abundance is the UV field, since the O2 
molecules are easily photodissociated, and the C/0 abun- 
dance ratio in the gas phase. As soon as C/0 is larger 
than 0.7, the O2/CO ratio falls below 0.1. In the frame of 
steady-state chemistry in a dark cloud, our O2 limit im- 
plies a gas phase C/0 ratio larger than 1 (at the la level) 
or C/0 0.9 at the 3(t level. This would mean that most 
of the oxygen is frozen onto grains. Another possibility 



Table 2. Derived upper limits at Icr 







N(02) 


N(CO) 


N(02) 


O2/CO 






NRO 


IRAM 


KP 


Best 


i^(GH 


z) 


56.26476 


136.845411 


118.7503 




r(l<T) 




<.025 


1500 


<0.12 




= 


5K 


<1.75el6 


1.51el9 


<3.47el6 


<1.16e-3 


Tx = 


lOK 


<4.57el6 


2.41el9 


<1.08el7 


<1.89e-3 


Tx = 


15K 


<8.90el6 


3.91el9 


<2.24el7 


<2.20e-3 



is to relax the hypothesis of steady-state chemistry, and 

consider time-dependent models (Prasad & Huntress 1980, 
Graedel et al 1982, Langer et al 1984): much more atomic 
carbon and oxygen exists at times well before steady-state, 
at so-called 'early-times', i.e. 1-3 10^ yrs. If steady-state 
has no time to establish, because of a short dynamical 



mixing time between the diffuse photodissociated phase 
and the dense dark phase, the relative abundance of O2 
can be reduced by one or two orders of magnitude (Chieze 
ct Pincau dcs Forcts 1989). Most of the oxygen would then 
be in atomic form, as supported by the recent detection 
of [01] 63/U in absorption towards Sgr B2 (Baluteau et al 
1997). Also in the model of Le Bourlot ct al. (1995) molec- 
ular clouds reveal chemical bistability, and there may exist 
a high ionization phase of the ISM with a very low oxygen 
abundance, where the O2/CO ratio is lower by two orders 
of magnitude than in the low ionization phase. 

The interstellar abundance of molecular oxygen re- 
mains problematic; balloon and satellite experiments are 
much needed. If confirmed, the present low upper limit 
for the O2/CO ratio strongly constrains chemistry mod- 
els. If a large fraction of the oxygen abundance is frozen 
onto grains we expect enhanced O2 and H2O abundances 
in regions of shocked molecular gas. 
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